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We study the electronic properties ol electrons in flat and curved zigzag graphene ribbons using 
a tight-binding model within the Slater Koster approximation. We find that curvature dramatically 
enhances the action of spin orbit effects in graphene ribbons and has a strong effect on the spin 
orientation of the edge states: whereas spins are normal to the surface in the case of flat ribbons, 
this is no longer the case in the case of curved ribbons. We flnd that for the edge states, the spin 
density hes always in the plane perpendicular to the ribbon axis, and deviate strongly from the 
normal to the ribbon, even for very small curvature and the small spin orbit coupling of carbon. We 
find that curvature results also in an effective second neighbor hopping that modifies the electronic 
properties of zigzag graphene ribbons. We discuss the implications of our finding in the spin Hall 
phase of curved graphene Ribbons. 

PACS numbers: 



I. INTRODUCTION 

The electronic structure of graphene depends both on 
its structure at the atomic scale, determined by the sp^ 
hybridization, and on its structure at a much larger 
length-scale, determined the shape of the sample-^-. 
Thus, the electronic properties of flat graphene differ in 
subtle but important ways from rippled graphene'^ and 
the properties of carbon nanotubes are fully determined 
by the way they fold^. Curvature is believed to affect 
transport^i^, magnetic^ and spin relaxation properties of 
graphene^. In a wider context, the interplay between 
mechanical deformations and electronic properties, the 
so called flexoelectronics, is giving rise to a new branch 
in nanotechnology. Whereas conventional electronics de- 
vices are based on the capability to tune their working 
properties by application of external perturbations in the 
form of electric and magnetic fields, mechanical deforma- 
tion can have a major impact on the properties of nano- 
electronic devices. This results in a wide range of new 
effects, like piezoelectric nanogenerators- and field effect 
induced by piezoelectric effect^ in ZnO nanowires and 
stress driven Mott transition in VO2 nanowiresi^. 

The spintronic and magnetic properties of graphene 
are intriguing. From the theory side, there are two bold 
predictions. First, graphene should display a Quantum 
Spin Hall phase with spin-filtered states in the zigzag 
edgesiiii^. Second, the same edges should have fer- 
romagnetic order—"—. Whereas the striking progress 
in the fabrication of atomically precise graphene rib- 
bons has not permitted to confirm the presence of sta- 
ble zigzag edges ye t^'^i^^ , several groups have reported 
the observation of ferromagnetic order in graphene and 
graphite^^"— , in most instances in samples that contain 
many flakes, structural disorder or have been irradiated. 

The observation of the spin Hall phase in flat graphene 
would require reducing the temperature below the spin- 
orbit induced gap, which is smaller than 10/xeV— "— . 
From this point of view it would be desirable to increase 



the strength of spin orbit interaction in graphene. Hints 
of how this could be achieved come from experiments. 
On one side, the spin relaxation time of graphene, as 
measured in lateral spin valves, is in the range of 100 
ps^'^, much shorter than expected from the small size of 
spin orbit and hyperfine nuclear coupling^. Thus, some 
mechanism enhancing the strength of spin orbit interac- 
tion must be at play in these samples. A possible candi- 
date could be either curvatureii^^ induced by ripples and 
adatoms^i^. Curvature has been shown to enhance the 
effect of spin orbit coupling in the case of carbon nan- 
otubes, for which recent experimental work has reported 
zero field splittings induced by spin-orbit coupling split- 
tings in the range of 200 jieV-^ . 

The recently reported fabrication of curved graphene 
ribbons by unzipping carbon nanotubes^"— opens the 
way towards the experimental study of the effect of cur- 
vature on the edge states of graphene ribbons. Here we 
study this system from the theoretical point of view and 
we compare the spin properties of a graphene ribbon both 
for flat and curved ribbons. In flat graphene the tt bands 
are decoupled from the a bands, unless spin orbit cou- 
pling is considered. However, the effect of spin orbit cou- 
pling on TT bands occurs only via virtual transition to 
higher energy a bands. 

In the case of flat graphene, it has been verified that 
the effect of spin orbit on the tt bands can be properly 
described by an effective spin dependent second neigh- 
bor hopping between the tt orbitals. This is the so called 
Kane and Mele model^iii^, which predicts that graphene 
is a Quantum Spin Hall insulator with a spin and val- 
ley dependent gap and peculiar spin-filter zigzag edge 
states^^'^^. In the case of curved graphene, tt and a 
orbitals are coupled, and to the best of our knowledge 
the validity of the Kane-Mele model has not been tested. 
This is why we adopt a different strateg y^'^^i"^^" — and 
we use a four orbital tight-binding model, which includes 
both the TT orbitals as well as the s, Px and Py orbitals. 

The rest of this manuscript is organized as follows. In 
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TABLE I: On-site energies and Slater-Koster parameters in- 
volving the same atoms, carbon-carbon interaction and two 
different atoms, carbon-hydrogen interaction. All the values 
are in eV. 

section pl|) we describe the tight binding method used in 
our calculations and review some general results about 
the spin properties of the system . In section (jIII[) we 
present results for the electronic structure of flat zigzag 
graphene ribbons and compare with those of the Kane- 
Mele model. In section (jIVp we address the main point of 
this work, the electronic structure of edge states in curved 
graphene zigzag ribbons. In section (|V|) we discuss our 
results. 



II. FORMALISM 

In this section we briefly comment the two different 
tight-binding approximations used to calculate the elec- 
tronic structure and we provide some theory background. 

A. Slater Koster approximation 

In most of the calculations in this work we use a multi- 
orbital approach, taking into account the 4 valence or- 
bitals of the Carbon atom, s,x,?/, z similar to that used 
in previous wor k'^^i"^^" — . Thus, counting the spin, the 
single particle basis has 8 elements per carbon atom. In 
addition, we passivate the edge carbon atoms with a sin- 
gle Hydrogen atom for which a single s orbital, with the 
corresponding spin degeneracy, is included. The matrix 
elements of the Hamiltonian are computed according to 
the Slater Koster approach considering only first neigh- 
bor hoppings. For simplicity we approximate the overlap 
matrix as the unit matrix. We model both the carbon- 
carbon and carbon-hydrogen hoppings of graphene with 
a set of tight-binding parameters derived by Kaschner 
et al^ from comparison with density functional calcula- 
tions. We show these parameters in table (|l|. 

Spin orbit coupling is treated as an intra-atomic po- 
tential: 

Vso^XS-J2Li (1) 

where A is the spin orbit coupling parameter, S is the 
spin operator and Li is the orbital angular momentum 
operator acting upon the atomic orbitals of site I. The 
representation of this operator in the basis x, y, z is pro- 
vided in the appendix (fS| . Whereas there is no consensus 



regarding the value of the atomic spin orbit coupling in 
carbon, the values reported in recent work range between 
A = 4 and A = 8 meV— 1^>2^. In this work we always dis- 
cuss our results for values of A in that range and, when 
some physical insight is gained by so doing, for values of 
A much above the realistic range. 

B. One orbital tight-binding model 

The low energy physics of most graphene based nanos- 
tructures can be described with a tight binding model 
with a single orbital per atom, which can be taken as 
a ^ = 1 atomic orbital projected along the local normal 
direction to graphene surface, the so called tt orbitals. 
From the discussion above, it is apparent that the atomic 
spin orbit operator mixes orbital states in the same atom 
with different values of m. However, in some instances is 
still possible to describe the low energy sector of graphene 
with an effective Hamiltonian governed by the by tt or- 
bitals in which spin orbit gives rise to spin dependent 
hopping termaiiii^. Wc express the effective Hamiltonian 
using second quantization operators c| ^ that create one 
electron in the atomic site / with spin a: 

"Ho = X! Tla^Ja'cl^^Cj.a, (2) 
/, J,(t,(t' 

The Hamiltoniam matrix is the sum of four terms: 

+ itRU ■ (tv,^' X di) n'^]] (3) 

The elements of the matrix Nj^^j ( Nj'^j ) are equal to 
1 when / and J are first (second) neighbor and zero ev- 
erywhere else. Thus, the first two terms are the spin 
independent first and second neighobour hopping. The 
third term is the Kane-Mele spin orbit modeUiiiSi. It is 
a spin dependent second neighbor hopping between sites 
/ and J which have a common first neighbor C. The 
unit vector along the bond between sites / and C (C and 
J) is denoted by di (^2)- In the Kane-Mele spin-orbit 
model the spin dynamics is linked to the bond orienta- 
tion. Thus, in flat graphene and graphene and graphene 
ribbons, the bonds lie in a plane so that the Kane-Mele 
spin orbit conserves the spin along the normal to the 
plane. This is in contrast to the curved ribbons and 
nanotubes considered below, for which the bond vectors 
are not restricted to a plane and no component of the 
spin operator is conserved. The last term in equation 
^ is the so called Rashba Hamiltonian, a first-neighbor 
spin dependent hopping. In contrast to all the other 
terms, the Rashba explicitly breaks inversion symme- 
try so it is associated to an external electric field. In 
the rest of this paper we calculate the band-structure of 
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graphene based one dimensional structures using the 4 
orbital Slater Koster model and compare with the re- 
sults of the 1 orbital model defined by eq. ([2|) and ([3|), 
taking the Rashba term equal to zero. Whereas the 1 
orbital model gives results very similar to those of the 
flat ribbon, this is not so in the case of curved ribbons. 



(a) (b) 




FIG. 1: Color online, (a) Scheme of the procedure to gener- 
ate curved ribbons, (b) Section of a cruved ribbon, (c) Per- 
spective view of a curved graphene ribbon with edge atoms 
passivated with hydrogen. 



C. Some general results 

We study one dimensional ribbons formed repeating a 
crystal unit cell. Taking advantage of the crystal sym- 
metry, the Hamiltonian can be written as '^'Hn.mik) 
where the indexes n and m run over the single particle 
spin-orbitals of the unit cell. For the one dimensional 
structures considered below, the unit cell is shown in fig- 
ure 1. The period of the crystal is given by the graphene 
lattice parameter, a, k is given by a scalar k and the 
Brillouin zone can be chosen in the interval —n/a and 
w/a. The eigenstates of the crystal are labelled with the 
band index and their crystal momentum k. They are 
linear combination of the atom / with quantum numbers 
n = l,m: 

^.kiN,r) = e'"'"' C,k{I,n,s)Mr- ri)Xs (4) 

/,n,s 

where N is an integer that labels the unit cell, and — 
77) is the atomic orbital with orbital quantum numbers n 
(which encodes I and to ) in the atom / of the unit cell. 
The eigenstate of the spin operator along the z axis is 
denoted by Xs where s can take values ±i. Because of the 
spin orbit coupling it is important to specify the positions 
of the atoms with respect to the spin quantization axis. 

Due to time reversal symmetry, every state with energy 
e^{k) must have the same energy that its time reversal 
partner, e^'i—k) where v and v' label states related by 
time reversal symmetry. In systems with inversion sym- 
metry the bands satisfy ty{k) = ey{—k) so that, in the 
same k point, there are at least two degnerate states. In 
systems without inversion symmetry, like the curved rib- 
bons considered below, a twofold degeneracy at a given k 



point is not warranted. In this non-degenerate situation 
we can compute, without ambiguity, the spin density as- 
sociated to a given state with quantum numbers v, k as: 

(4fc(/))= ClM,n,s)Ck,Al,n,s')Ss,s' (5) 

i,n,s,s' 

where S^^a' are the Pauli spin 1/2 matrices. 

In the cases with inversion symmetry, like the flat rib- 
bon considered below, for a given k point there are at 
least two degenerate bands. Thus, any linear combina- 
tion of states of the degenerate pair, 5'^^ and ^i^'k is 
also eigenstate of the Hamiltonian and has different spin 
density. In these instances, we include a infinitesimally 
small magnetic field in the calculation which breaks the 
degeneracy and permits to attribute a given spin den- 
sity to a given state. When the calculated spin densities 
so obtained are independent on the orientation of the in- 
finitesimally small magnetic field, they can be considered 
intrinsic properties of the spin states. As we discuss be- 
low, this is the case of the spin filter states in flat spin 
ribbons, which point perpendicular to the ribbon and 
have a strong correlation between spin orientation, edge 
and velocity, as predicted by Kane and Mel o^^i^^ . 

In order to characterize the properties of a given state 
it will also be convenient to calculate their sublattice po- 
larization: 

«fc>- E ICfc,.^'^,^)!'^.^) (6) 

/,n,s 

where az{I) = +1 when I is a A site and (Tz{I) = — 1 
when / is a B site. 



III. FLAT GRAPHENE ZIGZAG RIBBONS 
A. Two dimensional graphene 

The spin properties and electronic structure of the fiat 
ribbons considered below can be related to those of the 
two dimensional graphene crystal. We briefiy recall the 
the spin-orbit physics of two dimensional graphene as de- 
scribed within the Slater Koster modeP". Within this 
approach, the electronic structure of two dimensional 
graphene is described by a 16 by 16 matrix, correspond- 
ing to the two atoms A and B of the unit cell^S. At zero 
spin orbit the 16 bands are two copies (one per spin), of 3 
bonding-antibonding pairs of a bands and one bonding- 
antibonding pair of the tt orbitals which compose the 
states of the bands at the Fermi energy and are decou- 
pled from the a bands. The Fermi surface is composed 
of two points, K and K' , where the gap between the two 
TT bands vanishes. In the neighborhood of both K and 
K' points the two tt bands are linear and the kp theory 
is formally identical to that of masless two dimensional 
Dirac electrons. Spin orbit couples the a and tt orbitals, 
producing anti-crossings away from the Fermi energy and 
opening a gap at the K and K' points. 
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Within the one orbital approach the Hamiltonian of 
graphene reads (see appendix ([B]) for details): 

H = F(fcKl + G(fc)s.a. (7) 

where the a operators act on the sublattice space, A and 
_B, 1 is the unit matrix in the spin space and = ±1 
labels the spin. Here F{k) is proportional to t and cor- 
responds to first neighbor hopping whereas G{k) is pro- 
portional to tKM corresponds to the Kane-Mele second 
neighbor spin orbit coupling. In the K and K' points the 
F function vanishes but the spin orbit term gives rise to 
a gap A50 = G{K). As we show in the appendix, the 
function g satisfies -G{K) = G{K') = 3^3^ KM ■ Thus, 
because of the spin orbit interaction the states at the 
Dirac points in graphene are described by the effective 
Hamiltoniar^ii: 



^SOTzCTzSz 



(8) 



where = ±1 labels the valley quantum number, 
^so = 3^/3tKM- Thus, for a given spin orientation 
Sz — ±1, spin orbit opens a gap that takes a valley de- 
pendent value 2tzSz^so = i2Aso. This spin and val- 
ley dependent gap would make graphene a peculiar type 
of insulator which could not be connected with a stan- 
dard insulator by smooth variation of a parameter in the 
Hamiltoniar^iiii^, i.e., a Quantum Spin Hall Insulator. 



Flat ribbons 



with those at the Dirac points, that are described by the 
effective Hamiltonian Thus, as we move from valley 
K to K' (positive velocity bands) the spin t (|) states 
of the edge with sublattice A (B) must change from the 
top of the valence band to the bottom of the conduction 
band. The roles of spin and sublattice are reversed when 
considering the two bands that start at K' and end at 
K. Thus spin f electrons move with positive velocity in 
one edge and negative velocity in the other. 
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FIG. 2: (a) Bands for flat ribbon without SO. Inset: Zoom 
of the edge bands for case (a), (b) Edge bands with SO (A — 
SOOmeV) and Ikm = 0.42meV. (c) Slope vs A'^ (see equation 
). (d) Slope vs tKM 



We now discuss the electronic structure of fiat 
graphene ribbons with zigzag edges. The unit cell that 
defines the zigzag ribbon has N carbon atoms and 2 hy- 
drogen atoms that passivate the two carbon atoms in 
the edge of the ribbon. These structures were proposed 
by Nakada et a^. Using the 1-orbital tight binding 
model, without SO coupling, they found that zigzag rib- 
bons have almost flat bands at the Fermi energy, local- 
ized at the edges. An important feature of zigzag edges 
is the fact that all the atoms belong to the same sub- 
lattice. Since the honeycomb lattice is bipartite, a semi- 
inflnite graphene plane with a zigzag termination must 
have zero energy edge states^ whose wavefunction de- 
cays exponentially in the bulk, with full sublattice polar- 
ization. In finite width zigzag ribbons, the exponential 
tails of the states of the two edges hybridize, resulting in 
a bonding-antibonding pair of weakly dispersing bandsi^i. 
The bands of the zigzag ribbon can be obtained by fold- 
ing from those of two dimensional graphene, either with 
real or imaginary transverse wave vector, and the lon- 
gitudinal wave vector varying along the line that joins 
the two valleys, K and K' . Thus, the valley number is 
preserved in zigzag ribbons. 

Spin orbit coupling, described with the 1 orbital model, 
has a dramatic effect on the (four) edge band a^^'^^ . The 
second-neighbor hopping makes the single edge band dis- 
persive, and overcomes the weak inter-edge hybridiza- 
tion. Interestingly the quantum numbers connect well 



This scenario is confirmed by the 4 orbital model and is 
expected based upon the fact that 2D graphene with spin 
orbit coupling is a Quantum Spin Hall insulator. In figure 
2 wc show the bands of a flat ribbon with = 20 carbon 
atoms. In figure 2a we show the bands calculated without 
spin orbit coupling. The calculation shows both the edge 
and confined tt bands as well as some a bands higher in 
energy. In the inset we zoom on the edge states to show 
that they are almost dispersionless except when k gets 
close to the Dirac point. In figure 2b we show the same 
edge states, calculated with SO coupling, both within 
the 4 orbital and the 1 orbital model. For this particular 
case we have taken A = 500 meV and Ikm — 0.42 meV. 
Figure ([SJ:) shows a slightly different slope for valence and 
conduction bands for the four orbital case. This electron- 
hole symmetry breaking can not be captured with the one 
orbital Kane-Mele model. 

It is apparent that the edge states acquire a linear dis- 
persion e — mk. The slope m of the edge states dis- 
persion increases linearly with the size of the gap at the 
Dirac points, which in turn scales quadratically with A 
(and linearly with Ikm)- We can use m to quantify the 
effect of spin-orbit coupling on the edge states. Figure 
([2])d shows that, for flat ribbons, we can fit m = 

with a ~ 24^. 

Since the fiat ribbons have inversion symmetry, the 
bands have a twofold degeneracy. In order to avoid nu- 
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merical spin mixing of the degenerate states we apply a 
tiny magnetic field (always less than 2 • 10~^T) to split 
the states. As long as the associated Zeeman splitting is 
negligible compared to the spin orbit coupling, the direc- 
tion of the field is irrelevant. By so doing, we can plot 
the spin density of the edge states without ambiguity. In 
figure ^ we show both the spin density of the four edge 
states with k = 0.99-^ and the square of the wave function 
for valence states, calculated with the 4-orbital model. In 
agreement with the 1-orbital Kane-Mele model, the spin 
densities are peaked in the edge, oriented perpendicular 
to the plane of the ribbon. We have repeated the calcu- 
lation rotating the plane of the ribbon and obtained the 
same result. From inspection of figures ^jp) and dS]), it is 
apparent that the valence bands correspond to the edge 
states and display the spin filter effect, i.e., in a given 
edge right goers and left goers have opposite spin. 
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FIG. 3: Properties of the edge state valence band, (a) Square 
of the Wave function for valence states with k — -1-0.99^. The 
sublattice polarization is apparent, (b) Spin density for states 
with k — -1-0.99-^. (c) Spin density of time reversal symmetric 
state with k — —0.99-^. The spin density of the 4 states is 
polarized perpendicular to the sample. 



The unit cell of the curved ribbons is obtained as frac- 
tion of a (n,n) nanotube, with radius R (see figure 1). 
For a given nanotube we can obtain a series of curved rib- 
bons with the same curvature and different widths, 
W, or different number of carbon atoms N. We can also 
study ribbons with the same N and different curvatures 
i?, using a parent nanotube with different n. Our curved 
ribbons are thus defined by W and R or, more precisely, 
by n and N. 



A. Energy bands 

The energy bands of curved ribbons is shown in figure 
dH) for a ribbon with iV = 20 and R = 4.1nm. There 
are three main differences with the flat ribbon. First, 
the edge states are dispersive even with A50 = 0, as seen 
in figure Qb. This effect can be reproduced, within the 
one orbital model, including an effective second-neighbor 
hopping {t' ~ — 1.5meV). This dispersion breaks electron 
hole symmetry and competes with the one induced by 
SO coupling, as seen in panels (|31:,d,e), with A50 =5, 50 
and 500 meV respectively. 




IV. CURVED GRAPHENE ZIGZAG RIBBONS 

The calculation of the previous section , using the 4 or- 
bital model, backs up the conclusions of the Kane-Mele 1 
orbital model for the spin filter effect in graphene . How- 
ever, the bandwidth of the edge states is less than 0.1 fieV 
for the accepted values^ of A = 5meV. Thus, the effect 
is very hard to observe in flat graphene ribbons. This 
leads us to consider ways to enhance the effect of spin 
orbit. For that matter, we calculate the edge states in a 
curved graphene zigzag ribbon, similar to those reported 
recently^^'^^''*^, using the 4-orbital model. In contrast to 
flat ribbons, the properties of the edge states of a curved 
ribbon can not be inferred from those of a parent two 
dimensional compound, because it is not possible to de- 
fine a two dimensional crystal with a finite unit cell and 
constant curvature. 



FIG. 4: Electronic structure of curved zigzag ribbon with 
R = 4.1nm. (a) Bands for curved ribbon without SO. (b-e) 
Zoom of the edge bands for A = (b), A =5 meV (c) , A =50 
meV (d) and A =500meV (e). 

Second, curvature enhances the effect of spin orbit cou- 
pling, as expected. In order to separate the effect of spin 
orbit from the effect induced by curvature, we define the 
differential bands as the energy bands of the curved rib- 
bon at finite Xso subtracting the bands without spin- 
orbit: 

?,(fc) = e,(fc)-e,(A:,Aso = 0) (9) 

In figure ([5]) a we plot the differential bands for the rib- 
bon with = 20 and R —4.1 nm. They are two doubly 
degenerate linear bands with opposite velocities in the 
Brillouin zone boundary. Thus, when the effect of curva- 
ture alone is substracted, the dispersion edge states look 
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pretty much like those the flat ribbons in the region close 
to the Brillouin zone boundary. Thus, we can also char- 
acterize them by the slope of the linear bands, m. In 
figure (Ob we plot that the slope of the edge bands as 
a function of A , obtained with the procedure just de- 
scribed, for the same ribbon discussed before ( = 20 , 
R =4.1nm). It is apparent that the slope m is no longer 
a quadratic function of , in contrast to the case of flat 
ribbons. Even more interesting, in figure ([5]:) we plot 
the slope m for a fixed value of A =5 meV, as a function 
of the curvature k — R^^ . We find a dramatic 100- fold 
increase at small k. This result is consistent with the 
effective kp Hamiltonian for carbon nanotubes^i^. 




FIG. 5: (a) Differential energy bands (see eq. [9]) for ribbons 
with A'' = 20, R = A.lnm and A = 5meV It is apparent that, 
for small k the relative energy bands are linear, with slope m. 
(b) Slope of the relative energy bands as a function of A, for 
R = 4.1nm. (c) Slope of the differential energy bands as a 
function of k = R~^, both in linear and logarithmic scale, for 
A = 5meV. A 100-fold enhancement of the slope occurs in a 
very narrow range of small curvatures. 

The curvature induced enhancement of the spin-orbit 
effect on the edge states of the curved ribbon can be un- 
derstood as follows. In flat graphene the Dirac bands 
are linear combination of atomic tt orbitals with quan- 
tum numbers Z = 1, to = 0. The effect of atomic spin 
orbit coupling, XL ■ S can be understood perturbatively. 
To first order, spin orbit coupling has no effect on the 
product states |cr) x |/ = 1,to = 0) states. Second order 
coupling, via intermediate states with A with respect to 
the Dirac point, and orbital quantum numbers 1 = 1, 
TO = ±1 results in an effective spin orbit Hamiltonian 
acting upon the tt orbitals, with strength -j-, that con- 
serves Sz. Curvature changes this situation, because it 
mixes the tt orbitals with the Z = 1,to ^ orbitals, re- 
sulting in an spin orbit Hamiltonian for the electrons at 
the Fermi energjjii^ linear in the spin orbit coupling A . 

The third difference with the flat ribbon is apparent 
for the bands away from the zone boundary: they are 
not degenerate. This is shown in flgure ©a, which is 



a zoom of flgure Qb, including states at both sides of 
the Brillouin zone boundary. The lack of degeneracy is 
originated by the lack of the inversion symmetry of the 
curved ribbon. Interestingly, the degree of sublattice po- 
larization, (cTz), shown in flgureQb, anticorrelates with 
the splitting. In other words, the states strongly localized 
at the edges are insensitive to the lack of the inversion of 
the structure, which is non-local property. 
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FIG. 6: (Color Online) Detail of the edge states for curved 
ribbon with R = 4.1nm and A = 5meV. (a) Dispersion of the 
edge states (b) Sub-lattice polarization as a function of k 
for the two lowest energy bands. 



B. Electronic properties 

We now discuss the spin properties of the edge states 
in the curved ribbons. Notice that, due to the lack of in- 
version symmetry, there is no degeneracy at a given k so 
that the spin density is an intrinsic property of the state. 
In figure ([7]) we plot the magnetization density (S'i/,fc(/)) 
for the two lowest energy edge bands with momentum 
k = ^ + 0.01 (upper panels) and k = — 0.01 (lower 
panel), for a ribbon with R = A.lnm and A = bmeV. 
Whereas the correlation between the velocity, the spin 
orientation and the edge is the same than in flat ribbons, 
it is apparent that the quantization axis is no longer par- 
allel to the local normal direction. The spin of the edge 
states lies almost perpendicular to the normal direction. 
Thus, this is different from the case of nanotubes, where 
the spin quantization axis is parallel to the tube main 
axis, and different to the flat ribbon. 

The effect is even more striking in the case of an al- 
most flat ribbon, shown in figure ([5]) for which the spin 
quantization direction is clearly not perpendicular to the 
ribbon plane. Thus, a very small curvature is enough to 
change the spin quantization direction of the edge states. 
This is better seen in ^ , where we plot the angle formed 
between the spin quantization axis and the local normal 
in the edge atom. In (|ni)a we show the evolution of the 
angle for the states with k = —- -1-0.01 for the two lowest 
energy edge bands. For flat ribbons R~^ = 0, the angles 
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are and 180, i.e., the quantization axis is pependicu- 
lar to the ribbon plane. In the opposite hmit, for large 
curvature, the spin quantization angle lies perpendicular 
to the normal, i.e., tangential, but always in the plane 
perpendicular to the ribbon transport direction. 




FIG. 7: Spin densities for states with k — — +0.01 (panels 
(a) and (b)) and it = ^ - 0.01 (panels (c) and (d)). 




as a function of curvature k. For small curvatures the 
derivative blows up exponentially. Thus, the spin orien- 
tation of the edge states is very sensitive to moderate 
buckling deformation of the ribbon. This sensitivity is 
specially important for the small values of A adequate 
for carbon. Larger values of A reduce the effect. We have 
verified that the single orbital model with the general- 
ized Kane-Mele Hamiltonian is not sufficient to capture 
the effect. 




R '(nm 

°0 0.5 1 1.5 2 °0 0.05 0.1 

1/R(nm"^) 1/R(nm"^) 

FIG. 9: Angle 9 formed by the spin density and the local 
normal as a function of curvature k = for A = 5 meV. 
Notice the extremely large 



V. DISCUSSION AND CONCLUSIONS 




FIG. 8: Spin densities for states with k — + 0.01 (panels 
(a) and (b)) and k = ^ — 0.01 (panels (c) and (d)) in almost 
fiat ribbon, = 3.7- 10~'^nm~^. Notice that curvature can 
not be appreciated. 

The transition between the two limits is far from 
smooth. Even for the less curved ribbon that we have 
considered, with R — 271nm, we have 6 ~ 70°. This 
is better seen in where we show the low curvature 
region only for one of the bands. The dramatic effect 
of curvature on the spin orientation of the edge states is 
quantified in the inset of figure where we show 



We are now in position to discuss wether or not the 
spin filter effect could be more easily observed in curved 
graphene ribbons. One one side, the bandwidth of the 
edge states is dramatically increased. For moderate cur- 
vatures (i? = 4.1nm) the bandwidth of the edge states 
is 20/ieV, and it can reach 60 neV for R ~ Inm, to be 
compared with 0.1 fj,eV for flat ribbons. Thus, experi- 
ments done at 100 mK could resolve the edge band of 
curved ribbons. On the other hand, the curvature in- 
duced second neighbor hopping breaks the electron-hole 
symmetry (see figure dS^)), so that now there are 4 edge 
bands at the Fermi energy. At each edge there would be 
left goers and right goers with the same spin orientation, 
although slightly different k. Taken at face value, this 
would imply that there is no spin current in the ground 
state and backscattering would not be protected. On 
the other hand, the lack of eh symmetry implies that 
the edges are charged. The inclusion of electron-electron 
repulsion, even at the most elemental degree of approxi- 
mation, will probably modify the bands and restore the 
electron-hole symmetry. This is the case in most DFT 
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calculations of zigzag ribbons. Such a calculation is out 
of the scope of this manuscript. 

In conclusion, we have studied the interplay of spin or- 
bit coupling and curvature in the edge states of graphene 
zigzag ribbons. Our main conclusions are: 

1. In the case of flat graphene ribbon, the microscopic 
4-orbital model yields results identical to those of 
the Kane-Mele 1 orbital model. In particular, the 
edge states have the spin filter propertjsiiii^ and 
the spin is quantized perpendicular to the sample 

2. Curved graphene ribbons also have spin-filted edge 
states states. The bandwidth of the edge bands 
of curved ribbons is increased by as much as 100 
for moderate curvatures and is proportional to the 
curvature, for fixed spin orbit coupling. 

3. Curvature induces a second neighbor hopping 
which modifies the dispersion of the edge states 
and, in this sense, competes with their spin-orbit 
induced dispersion. 

4. The spin of the edge states not quantized along 
the direction normal to the ribbon. For moderate 
curvature, their quantization direction is a func- 
tion that depends very strongly on the curvature of 
the ribbon. Above a certain curvature, the quan- 
tization direction is independent on curvature and 
perpendicular to both the normal and the ribbon 
direction. The strong sensitivity of the spin orien- 
tation of edge states on the curvature suggest that 
flexural phonons can be a very efficient mechanism 
for spin relaxation in graphene. 

Note Added: During the final stages of this work a 
related preprint has been posted^ 
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Appendix A: Atomic orbital basis and matrix 
elements 

In this appendix we give the expressions for the atomic 
orbitals and the corresponding angular momentum ma- 
trix elements, necessary to compute the spin orbit ma- 
trix. In our formalism, the atomic orbitals are described 
in terms of the cartesian basis, s,Px,Py,Pz, which is re- 
lated to the basis of eigenstates of and L^, through 



the following transformationii^ 

\s) = \l = 0,m = 0) 
1 



(Al) 



= --=(|/ = l,m = l)-|? = l,m = -l)(A2) 
y2 

\py) = -^(|/ = l,m=l) + |; = l,m = -l) (A3) 
\Pz) = |; = l,m = 0) (A4) 
The spin orbit Hamiltonian operator reads: 



V50 = A 



(A5) 



which only affects the I = 1 subspace. The matrix ele- 
ments of this operator in the cartesian basis read: 





\P.A) 


\PyA) 


\pzA) 


\Px,i) 


\Py,i) 


\Pz,i) 







-iX/2 











A/2 


\PvA) 


iX/2 














-iA/2 


|p.,t) 











-A/2 


iX/2 


(A6) 


\PxA) 








-A/2 





iX/2 





\Pyii) 








-iX/2 


-iX/2 








\Pz,V) 


A/2 


-iX/2 















Appendix B: Graphene in the 1 orbtal model 

In this appendix we describe two dimensional graphene 
within the 1-orbital spin orbit model. The two atom unit 
cell is taken dimer forming degrees with the horizontal. 
The atoms are labelled A and B. We use the crystal 
vectors ai = a (1,0)) and 02 = a (Cos-|, S'm-|) . The 
first neighbor hopping between the R atom of a unit cell 
and its first neighbours reads 



F{k) = < (1 



(Bl) 



which accounts for the fact that 1 of the firt neighbours 
is in the same cell. Two Dirac points at which F{K) — 
F{K') = vanishes are given by Ka — 0) and K'a = 
-Ka. 

For the Kane-Mele second neighbor spin-dependent 
hopping matrix elements involves coupling to 6 atoms, 
four of them in first neighbor cell and 2 of them in sec- 
ond neighbor cells. We have: 

G{k) = itKM (e**^''^' - e"^'^^ ^ gife-(S2-Si) _ ^^^^^ (22) 
It can be readily verified that G{K'a) ~ —G{Ka) = 

3^/3tKM- 
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